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Several phenomena related to the critical behaviour of 
non-interacting electrons in a disordered 2d tight-binding 
system with a magnetic field are studied. Localization 
lengths, critical exponents and density of states are com- 
puted using transfer matrix techniques. Scaling functions of 
isotropic systems are recovered once the dimension of the sys- 
tem in each direction is chosen proportional to the localization 
length. It is also found that the critical point is independent 
of the propagation direction, and that the critical exponents 
for the localization length for both propagating directions are 
equal to that of the isotropic system, v ~ 7/3. We also cal- 
culate the critical value A c of the scaling function for both 
the isotropic and the anisotropic system. It is found that 
AJf = vA? • Ac. Detailed numerical studies of the density of 
states n(E) for the isotropic system reveals that for an appre- 
ciable amount of disorder the critical energy is off the band 
center. 



I. INTRODUCTION 



The problem of Anderson localizationEJ m anisotropic 
sysytems has attracted considerable attentionaO re- 
cently. It is generally accepteda that anisotropy does 
not change the universality class and that the isotropic 
results are recovered once a proper scaling of the 
anisotropic results is performed. If the dimension of 
the system size is chosen to be directly proportional to 
the localization length, the system should be effectively 
isotropic. The difficulty in implementing such a proce- 
dure lies in the fact that the localization lengths are usu- 
ally not known a priori. It was found through detailed 
numerical calculations!] that this scaling indeed works. 
It was also shownQ that the probability distributions of 
the conductance in the two directions are exactly equal 
to each other, provided that the ratio of the sides of the 
rectangle is proportional to the ratio of the localization 
lengths in the two directions. These scaling results were 
obtained for an anisotropic system where all the states 
were localized. n 

It is well knownEl that non-interacting electrons are lo- 
calized in 2d disordered systems. There are, however, 
some exceptions to this rule. These include electrons 
having strong spin-orbit couplings, integer quantum Hall 
systemsEl, arid tight binding models with random mag- 
netic fieldso. The best known example is the integer 
quantum Hall plateau transition occuring in a 2d non- 
interacting system in a strong magnetic field. Extended 
states do not exist as a result of Anderson localization 
except at a singular energy near the center of each of the 



Landau subbandsH£9. The localization length diverges 
at these critical energies E c as £ cx \E — E c \~ v . 

Another important point is the universality of the con- 
ductance at the critical point jpf the Anderson transi- 
tion for the anisotropic systemcl. From the generalized 
scaling functions, it has been established that the ge- 
ometric mean of the the critical value A c of the scal- 
ing function A = (\m/M) (as a function of £/M) is a 
constant independent of the strength of the anisotropy. 
(Here Am denotes the finite size localization length of 
a quasi-one-dimensional strip of finite width M.) Nu- 
merical calculations in both twoH and threea dimensional 
disordered anisotropic systems support this claim. How- 
ever, the same is not true for the conductance. Numerical 
calculationso in three dimensional anisotropic systems do 
not support a universal value of the conductance for the 
geometric mean. This might be due to too small sizes 
used in the 3d system or to a lack of universality of the 
conductance at the critical point. 

In this paper we investigate the scaling properties of 
the finite size localization length Am and the critical value 
A c of the scaling function in a two-dimensional system 
described by a tight binding model in the presence of a 
magnetic field. Both the isotropic case, as well as the 
anisotropic case will be examined. This is perhaps the 
simplest system that exhibits the correct behaviour of 
the metal to insulator transition. To our knowledge, no 
such calculations have been previously reported for the 
anisotropic tight binding model with a constant magnetic 
field. Some of the questions we try to answer are: how 
does the anisotropy effect the critical behaviour, espe- 
cially, will there be one or two critical exponents for the 
localization lengths; how do the anisotropic quantities re- 
late to the corresponding isotropic ones, especially, can 
we expect the geometric mean of the two anisotropic val- 
ues to equal the isotropic value; and what are the values 
for the scaling functions at the critical point? In section 
H we describe the model and the numerical methods we 
used, in section ^ we present and discuss our numerical 
results and in section IV we summarize the conclusions 
of this work. 



II. MODEL AND METHODS 

In the tight binding model one has the Hamiltonian 



i 



where the summations run over lattice sites i and j. We 
consider only nearest neighbour interaction in the hop- 
ping integrals Vij. The effects of an external magnetic 
field, characterised by a vector potential A(VxA = B), 
enter the model via phases of the hopping integrals with 



the integral connecting lattice sites i and j by a straight 
line. In two dimensions with a magnetic induction B 
perpendicular to the plane of the system, one can choose 
the gauge for the vector potential in such a manner that 
the phases vanish in one direction within the plane and 
are integer multiples of some number 2na in the other 
direction, such that the value of a completely charac- 
terises the influences of the magnetic field on the system. 
In particular, the denominator of a rational a equals the 
number of bands in the density of states of the system 
without disorder. Introducing anisotropy into the system 
by choosing different amplitudes Uj in the two directions 
within the plane will effect only the position of these 
bands, not their number. We bring disorder into the 
system by independently choosing all the site energies 
from a rectangular distribution of width W centered at 0; 
thus W is a measure of disorder strength. Both W and 
E are measured in units of the largest hopping matrix 
element t, which is taken to be unity. 

As Ojur main method we use the transfer matrix 
methodcl, where a matrix Tm connects the amplitudes 
of a state at both ends of a quasi-one-dimensional strip 
of width M and length N ^> M. Due to the anisotropy, 
we have to do this in the two spatial directions seper- 
ately. Therefore we get two sets of parameters Xm,x and 
\m,v which lead to two seperate localization lengths £ x 
and ^ y in the x- and y-direction respectively. Scaling 
of the data is used to improve on the values of £ x and 
£ y , which are then analysed to find the critical energies, 
where the localization lengths diverge as well as the crit- 
ical exponents of these divergences. 

We obtain the density of states by using a Lanczos 
procedureEH to diagonalize the Hamiltonian on squares of 
(linear) size L. The energy level separation distribution 
function P(s) should deviate markedly from a Poisson 
distribution for the local density of states around the crit- 
ical energy, approaching the Wigner distribution for the 
unitary ensembletj'Ej. We also use the density of states 
to show that for sufficiently strong disorder the critical 
energy does not necessarily coincide with the band cen- 
ter. 



III. RESULTS 

To obtain the critical energy E c for the anisotropic 
tight binding model, first we calculate Xm,x and \M,y for 
different strip widths M and energies E above and below 
the critical energy E c . As the exact position of E c varies 



with the disorder strength W, the hopping integral t x < 1 
in the difficult hopping direction as well as the magnetic 
field parameter a, we restrict our investigation to one 
set of these parameters W — 0.1, t x — 0.8 and a = \. 
The data for the more localised states show that M/Xm 
versus M is a straight line. The inverse slope of each 
of these lines gives a first estimate for the localization 
lengths £ x or £ y respectively, thus the smaller the slope, 
the more extended are the corresponding eigenstates of 
the system. For energies closer to E c the lines would 
be essentially horizontal. In order to accurately obtain 
E c , we have systematically calculated Xm,x and Am, 9 for 
large M. The results are shown in Fig. [I], where we plot 
Ajvf,x and Xu,y for the anisotropic case for energies very 
close to E c . From Fig. [j] we can confirm the existance 
of an extended state. Notice that Xm/M decreases as 
a function of M, which signifies localized states. For 
localized states, Am eventually reaches its large-M limit, 
which is a constant, and therefore Xm/M decreases as M 
increases. However, as can be seen from Fig. [j], at the 
critical energy E c , Xm/M saturates to a constant due 
to the absence of length scales. For the case studied 
(W = 0.1, t y = 1.0, t x = 0.8 and a = 1/8) we find that 
the critical energy E c is between —2.966 and —2.965, 
but closer to the second value. From Fig. [I] we can also 
obtain the critical values A c of Xm /M for both directions 
of propagation. We find A* = 0.92 ± 0.01 and A^ 
1 ..')') ± 0.01 with a geometric mean of 1.13 ± 0.01. 

To confirm that the geometric mean A c of the two 
anisotropic values is related to the value for the isotropic 
case, we have calculated systematically Xm/M versus M 
for the isotropic system (W = 4.0, t x = t y = 1.0 and 
a = 1/8) for very large values of M. These results 
are shown in Fig. |^. From Fig. || we obtain that in- 
deed E c = T-3 . .d fl in this case, in agreement with pre- 
vious resultsli!3li-a that used different techniques to get 
E c . In addition, Fig. ||| shows clearly that at the critical 
point of the isotropic system A' c so = 1.10 ± 0.03, which 
is approximately equal to the geometric mean of the two 
anisotropic values A^ and h v c . 

The critical value of Xm/M is related to the expo- 
nent aa. that can be obtained from the multifractal 
analysisa of the cigenfunctions at the critical energy by 
A^ 1 = 7r(ao — d) where jL is the Euclidian dimension of 
the system. Huckesteirfj calculated A c ==-1.14 ± 0.02 
for a real space model, while Lee et al.EZI determine 
A c = In - (1 + \/2) ~ 1-13 for a network model, both of 
which are close to the value AJf° = 1.10 ± 0.03 obtained 
for the isotropic case of the 2d tight binding model with 
a constant magnetic field. 

The next step is to use the values for the localisation 
lengths obtained in this manner to plot Xm.x/M as a 
function of £, X /M and Xm^/M as a function of i y /M. 
After combining the data for all energies into one graph, 
one usually has to adjust the values for the localization 
lengths slightly to make the data fall on a smooth curve. 
Fig. U shows that these two functions are independent of 
the value of E, as expected for one-parameter scaling. 
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However, the two scaling functions differ in their large-£ 
limit: the value is higher for the easy-hopping direction. 
To compensate for this anisptropy effect we use the fol- 
lowing straightforward ideao: Xm,x (^M,y) is a length in 
the x— (y-) direction along the length of the strip, so 
the appropriate scale should be £ x However, M is a 
length measuring the width of the strip and therefore has 
to be scaled with the other localization length. Thus we 
plot ^f- • |f vs % and • |j vs in Fig. g Not only 
do we obtain the same scaling function for both, but it is 
also the same as the isotropic one which we included for 
reference. The isotropic case was for W = 4 and a = 1/8. 
Of course under the assumption of one parameter scaling 
the form of the (isotropic) scaling function should not 
depend on the values of W and a directly (as long as 
neither vanishes completely) but only parametrically via 
the loclaization length £(E, W,a). Thus, the product of 
the two rescaled anisotropic functions equals the square 
of the isotropic scaling function. Immediately it is seen 
from this that the isotropic scaling function equals the 
geometric mean of the two anisotropic scaling functions 



<V/\£ 



AI.x 



A 



M,y 



M 



M 



(3) 



as the rescaling factors £ x and £ y cancel each other. As 
we have shown before in Fig. |l| and Fig. 0, indeed Equ. 
U is obeyed. 

The procedure of fitting the data to a smooth scaling 
function provides us with more accurate estimates of the 
localization lengths which we can now use to determine 
the critical behaviour of £. In Fig. || we plot the local- 
ization lengths as a function of energy. One can clearly 
see that the states are less localized in the easy hopping 
direction, as was to be expected. Fig. ^ also allows us to 
estimate E c , the energy where the localization length di- 
verges. We expect this critical energy to be independent 
of the strip orientation, as a higher-dimensional system 
would undergo a phase transition at this point, and our 
data give a strong indication that E c is indeed the same 
for both directions. We estimate E c w —2.965 ± 0.001. 
This is consistent with the results obtained in Fig. [l| 

The divergence of the localization length near the crit- 
ical energy is expected to follow a power law 



=Ho\E-E c 



(4) 



with some critical exponent v. To test this hypothesis, 
we plot the logarithm of £ vs the logarithm of \E — E c \. 
The result is shown in the inset of Fig. ||. That our data 
follows a straight line rather reasonably reconfirms our 
estimate for E c , as the plot obviously is quite sensitive 
to the choice for that value. Furthermore, both sets of 
data can be fitted by the same straight line, giving the 
same critical exponent v « 2.3 ±0.1. Once again, this 
is the same as the isotropic n value and very close to the 
theoretically predicted valud3 of 7/3 for the isotropic sys- 
tem. 



The distribution of energy level separations in a given 
energy interval depends on the typical extension of the 
eigenstates of the system with eigenvalues in that energy 
region. Spatial overlap of eigenfunctions close in energy 
helps delocalizing the particle. In a finite system, more 
of the strongly localized eigenfunctions can be accomo- 
dated without significant overlap. The more extended 
the eigenfunctions become, it gets more and more dif- 
ficult to fit several into the finite space and they must 
be seperated in energy. This leads to the phenomenon 
of level repulsion, known from chaos theory. The corre- 
sponding distribution of level separations Si = Ei — Ei-\ 
goes to zero for small s. In contrast, the distribution for a 
range of localized eigenstates has a maximum at vanish- 
ing level separation. More specifically, Random Matrix 
Theory predicts^ a Poisson distribution for the localized 
case, and a Wigner distribution for the extended case. 
We have calculated the distribution of energy level sepa- 
rations p(s) for the anisotropic system studied in Fig. g. 
We find that for an energy range close to E c = —2.965, 
p(s) is Wigner-like, whereas for the other energy ranges 
it is Poisson-like. 

Level statistics for the isotropic svsteaa have been ex- 
tensively studied by Potempa et alJiJ'Ej and Batsch et 
aLEaEZI, proving the validity of the approach in distin- 
guishing localized from extended states. In addition, 
the level number variance T, 2 ({N)) — \ ' (N) ias been 
numerically obtained for the isotropicj-systemES, using 
the Chalker-Coddington .network modelEj, and compared 
with analytical theoriesE3 which give for the spectral 
compressibilty \ — (d — D(2))/2d, where D(2) is the 
multifcactal exponent of the wavefunction at the critical 
pointH Klesse and Metzler obtain x = 0.124 ± 0.006E1 
Numerically obtained *alucs for D(2) include 1.43Jp0.03 
for a continuum modelE3, 1.56 for a network mojdew, and 
1.62 ± 0.02 and 1.71 for a tight binding modelE-3'Eij. Due 
to the limited size of our systems we were not able to pro- 
duce results for x for our anisotropic model. The number 
of energy eigenvalues sufficiently close to the critical point 
is not large enough to give good statistics for the number 
variance E 2 ((iV)). This point has to be adressed in the 
future. 

Finally, Fig. || shows the positions of E c for isotropic 
systems at W = 2 and W = 4 to be different from the 
band center. Although a Gaussian is not the correct form 
for the density of states it is usually a reasonable fit. For 
the stronger disorder, W = 4, the best approximation is 
achieved with a gaussian centered at E — —3.7 with a 
standard deviation of a = 0.4. Fig. || strongly indicates 
E c = —3.4 (arrow in top panel of Fig. ||). Similarly, 
for the lesser disorder, W — 2, a gaussian centered at 
E = —3.38 with a standard deviation of a = 0.21. A 
plot similar to that for the more strongly disordered case 
shows that E c — —3.32. However, for W < 1 the critical 
energy E c lies at the center of the Landau band. 
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IV. CONCLUSIONS 

In summary, we have performed detailed numerical 
study of the scaling properties of highly anisotropic sys- 
tems in 2d, with a metal to insulator transition. Scal- 
ing functions of the isotropic systems are recovered once 
the dimension of the anisotropic system is chosen to be 
proportional to the localization length. It is also found 
that the critical point is independent of the propagation 
direction and that the critical exponents for the local- 
ization length in both propagating directions are equal 
to that of the isotropic system. The critical value A c 
of the scaling function for both the isotropic and the 
anisotro pic cases has been calculated. It is obtained that 
Aj?° = ^JK X C ■ Af = 1.10 ± 0.03. Finally, density of states 
calculations revealed that the critical energy lies away 
from the center of the Landau band. 
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FIG. 2. The critical value of the isotropic scaling function. 
From the large-M data we estimate it to be 1.10 ± 0.03. 
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FIG. 1. The critical values of the scaling functionns can 
be obtained from the large-M limit of Xm/M at the critical 
energy. From the almost symmetrical behaviour of the val- 
ues for the non-critical energies at either side of the critical 
one, we assume that the value for E c is between —2.966 and 
— 2.965, but closer to the second one. Top: difficult hopping 
direction; bottom: easy hopping direction. 
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FIG. 3. The scaling functions for the difficult (open sym- 
bols) and easy (filled symbols) hopping directions \m/M as a 
function of £/M. The localization lengths have been adjusted 
to better fit the data to a smooth curve. Energies are —3.0, 
-2.99, -2.98, -2.97, -2.969, -2.968, -2.967, -2.964, -2.96, 
-2.95, -2.94 and -2.93. 
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FIG. 4. Plotting the rescaled scaling functions (cf. text) 
for the difficult (open symbols) and easy (filled symbols) 
hopping directions together with the scaling function for an 
isotropic system (crosses). 
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FIG. 5. The divergence of the localization lengths at the 
critical Energy; open symbols: difficult hopping direction, 
filled symbols: easy hopping direction. The values are taken 
after the adjustments made to obtain Fig. ^. Inset: To ex- 
tract the critical exponent of the localization lengths we plot 
£ vs \E — E c \ in a log-log plot. Both exponents are found 
to be 2.3 ± 0.1, roughly equal to the theoretical value for the 
isotropic system. 
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FIG. 6. The density of states for the lowest subband for 
a disorder strength of W = 4.0 (top) and W = 2.0 (bottom), 
indicating that the critical energy is off the band center. A 
fit to a gaussian distribution suggests that the band center 
in the W = 2.0 case is at E w —3.38, whereas the critical 
energy is E c ~ —3.32. 
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